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^ ; Low-energy excitations in spin-| antiferromagnetic(AF) Heisenberg spin ladders are 
studied by bosonization and gauge -theoretical description. The AF Heisenberg mod- 
ON '. els on ladders are described by spin-| fermions and the systems are reduced to rela- 
tivistic gauge theories of fermions by a mean-field type decoupling and linearization of 
fermion's dispersion relation near Fermi points. There gauge field is nothing but phase 
degrees of freedom of "mean field" on links. It is explicitly shown that zero-mode part 
of boson fields for the bosonization of fermions plays an essentially important role to 
obtain correct results. In the 2-leg case, the lowest-energy excitations are spin-triplet 
magnon and they are described by three massive Majorana fermions. On the other 
hand in the 3-leg system, spin excitations on the top and bottom chains are described 
by two massless boson fields. It is predicted that if coupling between spins on the 
top and bottom chains is introduced, a phase transition occures at some critical cou- 
pling. Above the critical coupling, the system acquires an energy gap and low-energy 
excitations are spin-triplet magnon and spin-singlet excitation. 
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1 Introduction 



In the last few years, quasi-one-dimensional strongly correlated electron systems are 
one of the most interesting topics in the condensed matter physics. These systems 
are closely related to the high-Tc cuprates and studies on them may shed light on 
mechanism of the high-Tc superconductivity. Recently, however, it has been rec- 
ognized that, besides relationship with the high-Tc cuprates, electron systems in 
quasi-one-dimensional space exhibit very interesting behaviors themselves which are 
to be investigated |J. 

In this paper, we shall study spin-| antiferromagnetic(AF) Heisenberg models on 
ladders. At present it is well-established that quantum spin-| systems have gapless ex- 
citations on odd-number-leg ladders whereas only gapfull excitations on even-number- 
leg ladders. This fact has been shown by numerical calculations and experiments, and 
understood by analytical studies through the nonlinear-sigma model approach^. 

Bosonization is also useful for (quasi-)one-dimensional electron systems. Among 
studies by the bosonization, Shelton, Nersesyan and Tsvelik showed that massive 
excitation in the 2-leg spin ladder is described by massive Majorana fermionsf^]. They 
started with two decoupled Hubbard chains with belief that low-energy excitations 
of the repulsive Hubbard chain are described by the AF Heisenberg chain, and then 
studied effects of the inter-chain interaction between spins by perturbation. Later 
it was verified that their theory explaines experiments very wellQ. However, we 
encounter difficulties in deriving their results by straightforward use of bosonization 
formulas. Especially we cannot obtain their expressions for the staggered part of the 
local-spin density (Eq.(Al8) in Ref.||) which plays an essentially important role to 
derive the effective Hamiltonian of inter-ladder interaction (Eqs.(7) and (8)). This 
is one of our motivations to study the AF Heisenberg models on ladders in more 
systematic and detailed manner. 

We shall directly study the Heisenberg models though most of theoretical stud- 
ies are given for the Hubbard models^]. To this end we employ a gauge-theoretical 
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description of the AF Heisenberg model which is obtained by a mean-field type de- 
coupling with keeping phase degrees of freedom of "mean fields" as dynamical vari- 
ables. Very recently similar investigation on 2-leg quantum spin ladder was given by 
Hosotani|J. However, we think that results of that investigation are not physically 
acceptable. This point will be discussed later on. 

This paper is organized as follows. In Sect. 2, we shall give the gauge-theoretical 
description of the AF Heisenberg model on 2-leg ladder. Spin operators are expressed 
by spin-i fermion bilinears and decoupling of the intra-chain spin-spin interaction 
is performed. We keep phase degrees of freedom of link auxiliary fields which are 
regarded as gauge fields. We intergarte out high-momentum modes of fermions and 
verify that Maxwell term of gauge fields appears. Low-momentum effective theory is 
nothiing but a Schwinger model of multi-flavor fermions. Then by using bosonization, 
we shall study effects of inter-chain interaction rather in detail. There zero modes of 
boson field for bosonization play an essential important role, and because of them we 
derive the results in Ref.||. 

In Sect. 3, we shall study the 3-leg quantum spin ladder. Especially we forcus on 
low-energy excitations on the top and bottom chains. We use the same method of 
bosonization in the 2-leg case but we have to introduce redundant boson fields in order 
to separate physical excitations. To recover the original 3-leg spin system, we impose 
local constraints. We study how these constraints influence low-energy excitations 
by integrating out massive modes and show that two boson fields remain massless 
and they describe low-energy excitations in the 3-leg spin ladder. We also consider 
effects of spin-spin interactions (both ferromagnetic and antiferromagnetic couplings) 
between spins on the top and bottom chains and show that the system acquires 
an energy gap above some critical coupling. This prediction should be checked by 
numerical calculation etc. 

Section 4 is devoted for discussion and conclusion. 
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2 AF Heisenberg model on 2-leg ladder 

Hamiltonian of the / c -leg spin-ladder system is given as follows, 

lc ^ lc-1 ^ 

Hl c = ^ \j ^ Si^i jr Si jr + J' ^ Si :r Si yr +i , (2-1) 
i r=l r=l 

—* 

where site on Z c Teg ladder is labeled as (i, r — 1, l c ). The operator Si >r is the spin 
operator with the magnitude |, and we express it in terms of fermions i\)^ r = (ipf r ) 
(or = 1,2 or U), 

Si,r = ^Ir^hn ( 2 - 2 ) 

with constraint 

EV#Vfr = l> ( 2 -3) 

a 

where a is the Pauli spin matrices. In the expression in terms of the fermions, the 
Hamiltonian is invariant under a local "gauge" transformation, 

-> e'^W- (2.4) 

Actually quantum spin systems can be described by gauge theories. This was first 
shown for the chain case in Ref . [[/fl and very recently for ladder case in Ref . . However 
results in Ref. do not agree with those in Ref.0, and then we shall revisit this 
system in this section. 

The following identity is easily verfied 

Si ■ 4- = ~ [ - w&j, ^M) + 1 - {4a - wM- - 1)' • 

We consider the case J/J'>1 and start with decoupled Heisenberg chains. For each 
chain we introduce auxiliary link field U^ r and Lagrangian of the spin-chain system 
is given as (in the imaginary time formalism), 

(2.5) 
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We can determine amplitude of the link field U = \Ui ir \ by the mean-field calculation. 
In this paper we consider chain and ladder of finite length N and impose periodic 
boundary condition for the spin variables S_N_ r = SV r . By straightforward calcu- 
lation we obtain U — - for sufficiently large N. We shall take into account effect of 
phase degrees of freedom of Ui >T and parameterize it as 

U ir = -e iaA ^ ){i \ (2.6) 

71 

where a is the lattice spacing and A^} = (Aq, A^) is gauge field for gauge transfor- 
mation (|2~^ ). 

For Ui tr = U dispersion relation of the fermions is that of tight-binding model and 
Fermi points are located at k = ±kp = ±j- as the present system is at half-filling by 
the constraint ( [2.3[ ). We are going to focus on low-momentum excitations and inte- 
grate out high-momentum modes of fermions. To this end we introduce a momentum 
cut off parameter A(<C and employ the approximation of linear dispersion rela- 
tion for the momentum region \k ± kp\ < A. For low-energy excitations near the 
Fermi points, we define right-moving and left-moving fermions {ipRipi) as in the usual 
way. By integrating over the high-momentum modes there appear renormalization of 
Fermi velocity etc. However the most important term for low-momentum excitations 
is the Maxwell term of the gauge fields Ay . Therefore we include this term. After 
rescaling the fermion fields, the following 2-dimensional U(l) gauge theory is obtained 
as a low-energy effective field theory of the quantum spin chain 

4lL = -4^ + E ^rV^r" + \ A t^ D, = d, + iA^, (2.7) 

e a=l,2 a 

where the coupling constant e 2 is a function of J, a and A, 



e 2 



7T 

T 



4 J 

3 a 
4JA 1 



f-oA m i 



dx- 



cosa; 



-i 



3 / V a ln(f )/' 
A, (2.8) 



and 7-matricies are explicitly given by the Pauli matrices as 7 = a±, 7 1 = 02 and 
therefore 7 5 = 03. 



The Schwinger model (|2.7| ) can be solved by bosonization technique ||. Bosoniza- 
tion formula is given as 

Aua(t,x) = ^c , r ^e ±i{9 ^ +2 ^ {t± ^ )/L} iVo[e ±iv/ ^^ ( *' x) ], 
00 

<t>rva(t,x) = ^(a^ n e- 2 ™^/ L + h.c), (L = Na, u — L,R) (2.9) 

n=l 

where N m is the normal-ordering operator with mass m and commutation relations 
are [q rV a,Vr'u'a'] = i5 Tr >5 UI// 5 aa >i [afv,m a r'v>,m\ = 8rr'$aa'$vv'$nm and Cs are functions 
of p's which make the fermion fields ip rva {t,x) satisfy the canonical anticommutation 
relations (CAR). 

In the rest of this section, we shall consider the two-leg ladder case, i.e., r = 1, 2. 
Hamiltonian in terms of bosons is easily obtained from ( |2.7| ) and ( |2.9| ). As in the 
usual cases, it decouples into charge and spin parts. It is also shown that the gauge 
fields Aff appear only in the form of the Wilson line intergal Q r ,w — J-l/2 dxA{ 
and its conjugate P r ,w- 

Ho = H chain! + H c h a in2 

— Hqc + Hqs 



+ £ J_Z d * l N ^(^c + *l 2 c + f^c) (2-10) 

H 0S = 7^{Qls + Qlas} 

+ s£>5*(?*- + *-) < 2 - n » 



where c = aJ/ir and 



~^{4>rRa + (prLa), 
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-^\~4>rRa + 4>rLa) 



rC 



rS 



1 



$aC = -^($1C±$2C), a = ± 

similarly for <3> aC -, a iy etc, 

Qam PLam PRam 

PLam + PRam 7U = S,C 

and neutrality condition, which is obtained by varying the action with respect to Aq 
is given as 

Q +c = L/a = N 

Q-c = 0. (2.12) 

From (|2.10 ), it is obvious that excitations in the charge sector acquire mass by the 
gauge interaction and do not include low-energy modes. 

It is convenient to introduce fermions for spin and charge degrees of freedom, 



VI 



D ±i{qamv+2nPamv(t±x)/L} jy^ |g±i V^W>amy (t ,£) j 



(a = ± m = S,C v — L,R) 



where 



'\mv 



± 



>2mu 



Q±Cu = ^{<?M + Qiui ± (92I/T + *m)} 
<l±Su = t;{<?m - qiui ± [02v] - q2u[)} 



(2.13) 



(2.14) 



(2.15) 



likewise for p and 



1 /2 



(2.16) 



In Appendix A, C"s are explicitly given, and with them it is verified that VWi/s 
satisfy the CAR. Then it is shown that spin-part Hamiltonian ( [2.11 ) is nothing but 
that of free massless Dirac fields ip±s = {'<P±sr' 1 P±sl)- 

Let us consider the inter-chain interaction. Continuum limit of the spin operators 
Si r are given in terms of low-energy fermions as follows, 

S i>r — ► lS r (x), S r (x) = J r R{x) + J rL {x) + ( 2 - 17 ) 
where smooth part of the SU(2) currents are given as 

Xu =■ ^ru^ru (2.18) 

and staggered part is 

n r (x) =: ipl L ^iprR ■ + h.c. (2.19) 
In terms of these operators, the inter-rung Hamiltonian is expressed as 

Hmng = J'l J dx J\(x) • J 2 (x) + Hi(x) • n 2 {x) . (2.20) 



From (|2.9| ) and ( |2.10|) , it is easily seen that the J-term in ( |2.20|) is marginal whereas 



the n-term is relevant. Therefore we shall consider effects of the n-term first. In term 
of the fermions ip a mu defined by ( [2.131 ), the relevant interaction term is written as 



rrrel 

rung 



J'l J dx fii(x) ■ n2(x) 

J'l { 

— J dx(^l CL ip +CR + ip j _ CL ip_ CR -h.c.^ 



and the Hamiltonian under consideration is given by 

H T 2 el = H oc + H os + H™ l ng . (2.22) 

By introducing auxilary fields or mean-field type calculation for the four-Fermi in- 
teraction Hru l n g, we can determine the mass of ipamv Action with auxiliary fields is 
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explicitly given as 

J'l r 

S[ip, a] = S [ip] + i— J dx a(i)\ CL i) +CR + ip ] _ CL *p- C R - h.c.) 
J'l { 

+i— J dx a* ($.sitf+SR ~ iPUl^-sr ~ ^-sl^-sr - h.c.) 
J'l r 

+ —Jdxaa', (2.23) 

where SoM i s action for the Hamiltonian Hqc + Hqs- Gap equations which determine 
values of cr and a' are given as 

(^ +cl ^ +C r + ^Icl^-or ~ h.c.) CT0 = ia' , (2.24) 

(ip+SL^+SR ~ IplsL^SR ~ ^SL^-SR ~ h - c ->^ = «^0, (2-25) 

where (• ■ -) ao etc. denote expectation value evaluated with putting a = a etc. in 
S[i/),a] in ( |2.23|) . Evaluation of the LHS on ( p.25|) is straightforward in terms of 



Majorana fermion representation of Hqs- Majorana fermions are defined 



as 



C=^^-Su + ^- Sv ), P«=^:W+Sv-il>lsu), iy = R,L) (2.26) 



and 



= -i{Hpl P r)- E (ClCr)} 



n=l,2,3 



3 M sf , / Ac / . Ac xo 



where M s = ^a' . 

On the other hand for the charge sector, there already exists mass term of the 
boson fields $ a c as seen from ( |2.10| ). Therefore calculation of the LHS on ( |2.24| ) 
needs some approximation. Here we employ perturbative calculation in powers of o"o 



9 



assuming that the boson mass term $ 2 is more relevent than the cos term.Q In 

the study of the Schwinger model, vacuum expectation values (VEV's) similar to the 



above operators are evaluated ||10||. Here we use similar argument 



From the constraints fl2.12| ), the Hamiltonian Hqq is given as 
Hoc = e ^Pl W + ^{N* + (2 PR+c + N + ^Y} 



2 +W 1L- n 

TTC TTC . 0+L 

2l + 2l (2! "'- c + — 
Zero-mode part of the wave function of the above Hamiltonian is labeled by integers 



e 2 L o TTC TTC. @+VKn9 

H — tt P -w + 777 + 77t( 2 PR-c H — ) + nonzero modes. (2.28) 



(pr+c = n +,Pn~c = n~) and wave function is explicitly given as 

(q + ,q^\n + ,n^ = e -<fe+»++9-n-) _ (2 .29) 

If there is no mass term of fermions ip±cv, the above states |n + ,n_) are degenerate. 
The mass term, 

-iM c (^ Ci V+cfl + ^W-cfl-h.c), M c = ^a , (2.30) 
has nonvanishing matrix element between the states with different (n + ,n_): 

l^+c^+c-fll"^™-) oc<y J++ljn+ <5j_ )n _, (2.31) 
etc. Then we should define the ^-vacuum by taking linear combination of them, 

|0+,0_)= ^ e i0 +"+ +i0 - n -|n+,n_). (2.32) 

Including the contribution from the nonzero modes, we evaluate the matrix element 
in the 0-vacuum as follows, 



(^ +CL ^ +CR ) e = e-^e^J^ (2.33) 



1 This is certainly correct as seen from the conformal dimensions of these operators. For 
renormalization-group study of the massive sine-Gordon theory, see Ref . B . 
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etc., where 7 is the Euler constant. Then the VEV of the mass term ( |2.30| ) depends 
on the values of the 0's. It has the largest value for 9 + = 6- = ±7r/2, at which the 
ground state energy of the system has the lowest value. The final result is 



M -^f^ < 2 - 34 > 

Low-energy Hamiltonian of the spin sector is then given by 
H s = H s + H rung . s 

J -L/2 n =l,2,3 Z 

+{^(PLd x p L - PrOxPr) - 3iM s p L p R j\. (2.35) 

The same Hamiltonian with the above was obtained by Shelton et.al.|| from the 
Hubbard Hamiltonian, but as we mentioned in the introduction, straightforward use 
of the bosonization formulas ignoring the zero modes does not give (|2.35|) . From the 
derivation in this paper, it is obvious that the zero modes in the bosonization of 
fermions play an important role and we have to consider the ^-vacuum in order to 
have the correct ground state of the system. 

From (j2.35|) , low-energy excitations in the 2-leg quantum spin ladder are described 
by spin-triplet Majorana fermions £™ and spin singlet fermion p v . This was verified by 
experiments @]. Marginally relevant interactions like J\ ■ J 2 give only renormalization 
of velocity and mass M s . 

Recently Hosotani studied quantum spin ladder by a similar method with ours, 
but he did not obtain the above result 0. In his approach, nearest-neighbor (NN) 
fermions in the real space are paired in order to obtain relativistic Dirac fermions, 
though we have paired right-moving and left-moving fermions near the Fermi points 
in the momentum space, as it should be. We are not quite certain if real-space 
paring gives correct continuum limit which describes low-energy excitations. We also 
suspect that the bosonized expression of inter-chain interaction is not correct. This 
may be connected with the above problem. The inter-chain interaction should contain 
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not only the scalar fields (J) a but also the dual scalar fields <j) a because of the term 
(S1S2 + h.c.) in the Heisenberg Hamiltonian, but his expression does not. Because 
of this, his result does not seem to be SU(2) invariant. 



3 AF Heisenberg model on 3-leg ladder 

In the previous section, we studied low-energy excitations in 2-leg spin ladder. In 
this section, we shall apply similar analysis to the 3-leg case in which excitations are 
expected as gap less. Hamiltonian is given by (p7T| ) with l c — 3, and the same gauge- 
theoretical description is possible for each chain. Most relevent inter-chain interaction 
is described by a similar Hamiltonian to the n-term in (|2.20 ), 

H lfung = J ' 1 J dx (M x ) ' M x ) + M X ) ' M x )) ■ (3-1) 

Then it is convenient to introduce the following boson fields which describe excitations 
in the (1-2) chains, etc., 

¥±Au = {4>lu±4>2u)/V2, 

<P±Bv = (<f)3u±4>2u)/V2, (3.2) 

<P±Cv = (cj) 3u ±(j) lu )/V2, (v = R,L) 

where the zero modes are defined in the same way. We shall consider only spin 
excitations and have omitted subscript "S"', i.e., <f> r = (0 r j — <f> r ±)/y/2 (r = 1,2,3) in 
( ft .2D since charge excitations are massive as in the 2-leg case, 

i=l,2,3 X=A,B,C 

where notations are self-evident. It is also useful to introduce fermions *&±x (X — 
A, B, C) corresponding to the above boson fields for each channel. As we showed in 
the previous section, the fermions *&±a,b become massive as a result of the inter-chain 
interactions H^f ung . Spin part of the Hamiltonian of the 3-leg spin ladder is then given 
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by 

Hss= E (H +x + H- X ), (3.3) 

X=A,B,C 

where 

H +x = Jdx {i^ix L 9x^+XL - *\.xr9x*+xr) 

+ iM s (l - 6 Xt c)(*l XL $+xR ~ h.c.)} (3.4) 
H-x = Jdx{i~{*l XL d x *- XL -*l XR d x *- XR ) 

- iM s (l - 6x,c){*Ixl*-xr - h.c.) 

- 2*M S (1 - 6 x ,c)(*L XL *l XR - h.c.)}. (3.5) 

The A, B part of the above Hamiltonian is the same with that in the 2-leg case and 
M s = ±(J7/47r)e 7 v /2e 2 /7rc, where we set 6± A = e l ±6 2 = 9± B = 9s ±02 = ±f • The 
H± c is the system of massless fermions (i.e., H os ). In the above expression ( |3.3j ), ( |3.4j) 
and ( |3.5|) , we have neglected the marginally relevant interactions of fermions which 
appear from J-terms in the inter-chain interaction. These terms cannot be neglected 
in the present 3-leg system and we shall discuss it later on. 

From ( |3.2j ), it is obvious that not all the bosons tp±x or fermions *&± x are inde- 
pendent. We have to impose the following constraints: 

l P+Au — <f-Av = P+Bu — ¥>-Bv, 

¥+Bv + ¥-Bv = ¥+Cv + V-Cv, (3-6) 

ip + cu - v-cu = v+Au + (f-Au, [y = R, L) 

and similarly for zero modes. These constraints can be rewritten in terms of fermions 
by using the following identities, 

^±xl^±xl = P±XL/L + d x ip±xL/Vn, 

^±XR^±XR = P±XR/L + d x Lf ±XR /y/n. (3.7) 
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From dJ]) and fl3"77D,fl 



-Bv 



-Ci 



(3.8) 



The above constraints ( |3.8|) can be imposed by Lagrange multiplyers Aj„ (i = 1, 2, 3) 
and we add the following terms to the Hamiltonian H^f (|3.3|) , 



Through the above terms, A, B and C fermions are interacting with each other. 

We are interested in how the constraints (|3.8| ) influence low-energy excitations, 
i.e., the C fermions which are gapless if we ignore the constraints. To this end, 
we shall use a similar argument as in the study of the nonlinear-cr model, i.e., we 
integrate out the massive fermions ^±Av and ^±bv and see whether mass terms of 
the Lagrange multiplyers appear or not. If mass terms of Xi U are generated, the 
constraints ( |3.8| ) tend to irrelevant for low-energy excitations, whereas if \ v remains 
massless the constraint is intact. 

The action of the +A-sector S + a is written as 



2 Here we neglect the constraints on the g-terms in the bosonization. This is justified, since 
operator q is the conjugate to the momentum p and wave function is given as an eigen state of p as 
discussed in Sect. 2. 



= (Ai„ - \ Zv )¥ +Av ^ +A u ~ (Ai„ + \^ ] _ Au ^-Av 

+ (A 2 „ - \u,)*\. Bv *+Bw + (A 2 v + \lv)^-Bv^-Bv 
+ (A 3 , - \ 2v )^\ C u^+Cv ~ (A 3 „ + A 2l/ )tti c „tf_ C „. 



(3.9) 




(3.10) 
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where we put xq = cr/2 and 

A = -i(\ 1L - A 3 l + \ir - A 3jR )/2, 
Ai = —(Ail — A 3 l — Air + X 3R )/2. 

Then effective potential of the A^ is obtained as0 

-S$[A] = -i fd 2 x{^ +A {x)^A^ +A {x)) 



-- J d 2 xd 2 x (^ +A (x)^ fl A ll ^ +A {x)^ +A {x')-f II A fl ^ +A {x')) c 

/f d 2 p 

+ \j f d 2 x J j^Trlg+ip^A^g+ip^A^} 
J d 2 xa P (X 1L - A 3L - Air + \ 3R ) 2 (3.11) 



where 



9+ip) = {l^ + i{2M s /c)^)-\p) 



f + (2M 5 / C ; 



2 



(3.12) 



1 Ac/2M S 
ap = . (3.13) 

8^ yjl + (Ac/2M 5 ) 2 

The action of the — A-sector is given in terms of Majorana fermions defined by 

c 3 = + *i^)/2 

X, = (*-a„ - *l^)/2i, (3.14) 



and 



5-a = /rf 2 a:[-C 3 (7A + ^(2M s /c)75)C S 
1 

2' 

+ < 3 7 ^ X ], (3.15) 



+ T^xilA ~ «( 6M s/c)7s)x 



3 Plcase remark that regularization scheme is specified by the spatial momentum cutoff A. 
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where 



Bo = (Ail + A 3 l + \ir + X 3R )/2, 
Bi = -i(X 1L + X 3L -X 1R -X 3R )/2. 

Then contribution from the — A-section to the effective potential of the \ v is obtained 

as 



where 



S e Jl 



[A] 



-% \ d z x{C {x)l^ X {x)) 



~ J d 2 xd 2 x' (( 3 (x)^A„x(x)( 3 (x')^A,x(x')) c 



\!*>I ,P " 



2 

J d 2 x [b P {(X 1L + A 3L ) 2 + (Air + \ m ) 2 } 
+ c P {X 1L + A 3L )(Air + X 3R ) 



2 , /\ \ \2i 



(3.16) 



9i(p) = M + t(2M s /c) l5 y\p) 

p 2 + (2M s /c) 2 
gitp) = (7A-i(6M 5 /c) 75 )- 1 (p) 
iVvTlii -i(6M s /c)^ 5 



6 P = 



c P = 



p 2 + (6M s /c) 2 
1 Ac 



32tt 2M s 



{y/9 + (Ac/2M 5 ) 2 - ^/l + (Ac/2M S ) 2 ) 



3 , /Ac/2M 5 + Jl + (Ac/2M s f 
In - 



!6vr Va c /6M s + ^1 + (Ac/6M 5 ) 
Similarly contributions from +B, —£> -sectors are evaluated as 



qeff 



[A] 



J d 2 xa P (X 2 L - Ail - A 2 p + Aip) 2 , 
J d 2 x [b P {(X 2L + Ail) 2 + (A 2 p + X 1R ) 2 } 
+ c p (X 2 l + Ai L )(A 2 p + X 1R ) 



(3.17) 

(3.18) 
(3.19) 

(3.20) 



(3.21) 



(3.22) 
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Then the potential is given as 



qeff i qeff I qeff I qeff 
+A -+- D_ A -t- D +B -+- 



2(a P + b P )\\ + (a P + b P )X 2 2 + (a P + 6 P )A 



(2a P - cp)Air Ai - (a P - cp/2)A 2 r A 2 - (a P - cp/2)A 3 r A 3 



2(a P - 6 P )Ai(A 2 + A 3 ) + (2a P + c P )X 1 r (X 2 + A 3 ) 



(3.23) 



where A; = (X iR X iL ) and r = o x . 

As we mentioned before, there are also contributions from marginally relevent four- 
Fermi interactions which are generated by the smooth part of the spin interactions, 
renormalization effect of high-momentum modes, etc. We first ignore their effects. 
Then we intergate over the field Ai in S e ^ and obtain 

1 



where 



s' eff [x 2 ,x s 



d x 



^(X 2 + X 3 )(r]-ujT )(X 2 + X 3 ) 



+ -(A 2 -A 3 )(l-ar )(A 2 -A 3 ) 



(3.24) 



a = (ap - c P /2)/(a P + b P ) 

(3 = V2(a P -bp)/(a P + bp) 

S = V2(a P + cp/2)/(a P + bp) 

rj = \ + {f3 2 + 5 2 - 2a(35)/2{l-a 2 ) 

uj = a- {a({3 2 + 5 2 ) - 2[35}/2(l - a 2 ). 



(3.25) 



Then action of the C-sector is obtained as follows: 



S+c — 



S. 



c 



J d 2 x[ty +clil d^ +c + (A 2 „ - X 3u )¥ +Cu ^ +Cu 
-(A 2 -A 3 )(l-aT )(A 2 -A 3 )], 



d 2 x 



^{X 2 + X 3 ){r]-ujr )(X 2 + X 3 ) 



(3.26) 



(3.27) 



17 



By straightforward calculation, it is verfied that the mass matrix of the fields A2,3 in 
Q3.261 ) and (|3.27|) has positive eigen- values for the coefficients (|3.25| ), and then A2,3 



can be integrated out. Marginally relevant terms appear like ^/\_ cr ^ + cb.^+cl^+cl 
but the C-sector still remains gapless. The investigation given so far shows that spin 
excitations on the top and bottom chains are gapless and described by two massless 
boson fields <p±cv, though excitations in the middle chain have a gap. 

We shall comment on effect of the marginal interactions which we have ignored 
so far. As in the 2-leg case, the SU(2) current I v is given as follows in terms of the 
spin-triplet Majorana fermions, 

C = (*+a„ + *U) A C = (*+Av - *U,)/2i, (3.28) 

4, = ^ m "CC\ (3-29) 
where e lmn is the antisymmetric tensor. Then interaction terms are given as 

HT r = fdx [ 9l £ AbJal + 92( E CrCDXrXl] , (3.30) 

l n=l,2,3 

where g\ and gi are coupling constants. Similar interactions appear for the (2-3) 
chains i.e., the ±B sector and also terms like IapJbli e ^ c - ^ i s obvious that these 
interactions also generate mass matrix of the multiplyers through bubble dia- 
grams. We do not give detailed calculations of these diagrams but we believe that in 
wide region of the parameter space of #2) existence of the interactions (|3.30|) etc. 



does not change the above result and excitations on the top and bottom chains are 
described by the two massless scalar fields (p±c- 

Let us consider how the gapless modes on the top and bottom chains evolve in the 
presence of inter-chain interactions between these chains, i.e., we add the following 
term to the Hamiltonian, 

H'rung = J" X/ ^,1^,3- (3.31) 

i 

From the discussion given so far, one may think that the gapless modes acquire 
a gap and low-energy excitations are described by spin-triplet Majorana fermions. 



However here again careful study on the zero-modes of the bosons of the bosonization 
is required since the #-terms of the A, B, C-sectors are not independent with each 
other. It is easily proved that #'s must satisfy constraints similar to Eq. fl3.6D - 



9 + a — 9 -a = 9+b — 9-b 



9 + b + 9-b = 9 +c + 9_c, (3.32) 
9+c — 0-c — 6+a + 9-a- 



From ( |3.32| ), it is obvious that for 6±a = 9± B = ±| we have 9± c = (mod tt) and 



therefore the fermions $?±c remain massless in the ^-vacuum for J" = even if the 
interaction (|3.31| ) exists. 

It is expected that the ^-vacuum for J" = tends to unstable because of the 
interaction (|3.31|) and there exists some critical coupling J" at which values of the 
parameter #'s of the lowest energy change from those of the J" = case. In that 
new phase, fermions ^/±c acquire a mass. The critical coupling J" can be estimated 
from the low-energy effective action of the 3-leg system similar to Eq. fl2.23j ) of the 
2-leg case. We expect that the charge sectors are less sensitive to the J" coupling 
than the spin sectors since the charge-sector scalar fields are already massive and 
their leading-order contributions to the effective potential of cr's are cancelled by the 
quadratic terms of cr's similar to the last term in ( |2.23|) oo' . Therefore let us study 
effective potential obtained from the spin sectors and get a rough estimate of J". 

The effective potental of the auxiliary fields a' x (X = A, B, C) or the masses 
Mx are obtained by integrating over the fermions ^±x- To this end we ignore the 
constraints ( |3.8| ) as we can expect that these constraints are irrelevant for low-energy 
excitations. From (|3.32j ) independent variables are, e.g., 9 + a,9 + b and 9 +c . We 



futhermore require the symmetry between the first and the third chains, i.e., we take 
9 + a = 9 + b- Then the effective potential is obtained as follows by intergating over the 
fermions $±a,b,c] 



^ /KV8 . )=6( ^ r[ln( ^ r _ l]+2( ^ )1ln( 



V Ac 
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+ 3 (^7) I ln b7) + 



3lA( V2r ln( 3Mc f 1 



Ac 



Ac 



(3.33) 



where 



M AB = ^ • (sin 9 +A + sm(-8 +B + 9 A 



c, 



M C /M AB = (j"/f) 



sin 9 



+c 



(3.34) 



sin 9 +A + sin(-6» + B + 9 +c ) ' 

From ( |3.33j ), ( |3.34| ) and the fact that f(x) = x(lnx — 1) is a decreasing function for 
x < 1, it is obvious that V e ^ has the minimum at 9 +A = 9 +B — ±f and 9 + c = ±7r 
for J" = as in the 2-leg case. We shall study the stability of this ^-vacuum in the 
presence of the J" coupling. 

We have numerically studied the effective potential V e ^ in the whole region of 9 +A 
and 9 + c and found that the absolute minimum stays at (9 +A , 9 + c) = (f, 7r) for small 
values of J". However at some value of J", the absolute minimum starts to move 
away from that point continuously. In order to study this behavior more precisely, 
let us introduce variable <5's and parametrize #'s as, 



9 +A = rt/2 + 5 A , 9 +B = ti/2 + 6 B , 9 +c = 7i + 6 C , &a,b,c < 1. 



(3.35) 



We assume (A — B) symmetry, i.e., 5 A = 5 B . Then the masses are given as a function 
of <5's as 



Mab oc cos(5 a ) + cos(5 A - 5 C ) ~ 2 - (6 A - y) - 
Mc oc — sin(5c). 



(3.36) 



From ( |3.36| ), it is expected that the absolute minimum exists close to the line 5 A = 
5c/2. Then let us study V e ^ on the line 5 A = Sc/2. Numerical calculations are shown 
in Fig.l. 

The calculation of the effective potential in Fig.l indicates that transition to the 
phase with nonzero 9±c (mod 7r) actually occures at a certain valu of J". Though this 
result is obtained by a rough estimate of the effective potential as we have ignored 
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contribution from the charge sectors and the constraints ( |3.8| ), we expect that there 
really exists some critical value of J" at which phase transition occures and fermions 
ty±c become massive. It is possible that J" = or J" = J' for some values of J. If 
J" < J', low-energy excitations are spin-triplet magnon and spin-singlet excitation 
and ratio of their gaps is 3. 

Very recently symmetric antiferromagnetic case of the above system (i.e., J' = J") 



was studied by Kawano and Takahashi in Ref . [jTTJ and they obtain the result that the 
system has an excitation energy gap for any finite value of J' = J". This result 
supports the above prediction. 

Let us take a brief look at the strong-inter-chain coupling to observe the above 
transition in this limit. We assume J' , J" > J and concentrate on the states on each 
rung. There are three sites and totally eight states on each rung. Eigen states and 
eigen values are given as follows; 

J" 

\ot)x\S) 2 z - \S) 12 \a) 3 , E = -J' + —, 
3 J" 

\S) 13 \a} 2 , E = — — , (a =1,1) 

ITT!)i23 + IT!T)i23 + |ITT)i23, (T<-i), E = y + X' ( 3 - 37 ) 

J' J" 

I TTT)l23, | 111/123, E = — + — , 

where \S)ij = | l)j — \ f)j and the other notations are self-evident. Therefore 
for J' > J", the lowest-energy states are doubly degenerate and they have spin-| 
degrees of freedom. At finite J coupling, this spin-| degrees of freedom on each 
rung behave like spin on a single spin chain. This is the reason why the system has 
no energy gap (in the strong-inter-chain couplig limit) and behaves like the spin-| 
chain. On the other hand at J' = J" , degeneracy of the lowest-energy states doubles 
because of the permutation symmetry. Therefore at finite J, one can expect a phase 
transition near this point. Actually the system is maximally frustrated at J' = J" 
and the lowest-energy state has rather complex structure [JTTJ] . Our analysis in this 



paper suggests that this transtion remains in the weak-inter-chain coupling region. 
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4 Conclusion 



In this paper, we discussed low-energy excitations in the AF Heisenberg model on 2-leg 
and 3-leg ladders. We used gauge-theoretical description of quantum spin chain and 
bosonization method. It was shown that in the 2-leg system low-energy excitations 
are described by spin-triplet and spin-singlet Majorana fermions. In order to obtain 
this result, the zero- modes in the bosonization play an important role. Our study 
justifies the results obtained by Shelton et.al. in Ref.||. Then we apply the same 
method to the 3-leg case and showed that excitations on the top and bottom chains 
are described by two massless scalar fields. 

One may wonder how our method is applied to the 4-leg system. In order to obtain 
a low-energy effective field theory, we first integrate out high-momentum modes in 
the original lattice system. In the 2-leg and 3-leg cases, this procedure is rather 
straightforward. However in the 4-leg case, it is naturally expected that nontrivial 
renormalization of the inter-chain couplings appears, i.e., the couplings J\i and J34 
develop and J23 tends to small. In low-energy states, spin-singlet pairs form on the 
(1 — 2) and (3 — 4) ladders respectively. Then it is natural to describe low-energy 
excitations on the 4-leg ladder by starting with the system of two decoupled 2-leg 
ladders. At present this problem is under study and results will be reported in a 
future publication. 
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Appendix 



A Explicit form of C's 

In this appendix, we shall give one of the possible choice of operators C which appear 
the bosonization formulas ( |2.9| ) and (|2.13|) , 

Ci m = exp[ni(p 1L1 + p lL[ + p 21 + P214) + ^ip XL \\ 

= exp[2rrip^ C L + 2m(p +C L + P-cl + P+sl + P-sl)], 

= exp[iTi(pm + p m + pat + P2L1) + m{piRi + P2R1)} (A.l) 

= exp[27up_ CL + 7ii(p +CR + P+sr)), 

C 2 ri = exp[ni(pm + p 1Ll + p 2 \ + p 2 Li) + m(pm + Pm - P2L] - P2L1)} 

= exp [2nip_ CL + 2mp_ CR } , 

C\L\ = 1, 
C\Ll = C 2 L1 

= exp[ni(p 1Ll +p2Lt)] 

= exp[iri(p +CL + P+sl)], (A.2) 
C 2Ll = exp[ni(pm + p 1Ll - p 2L} - p 2Ll )} 

= exp[27rzp_cL]. 
Then for spin and charge fermions, 

C +CR = exp[m(p +CL + p_ CL + p +SL + P-sl) + i^P-CL 

+ ni(p +C R + P-cr + P+sr)}, 
C- C R = exp [ni (p +C L + P-cl + P+sl + P-sl) - ^P-cr), 
C+sr = exp[iri(p +C L + P-cl + P+sl + P-sl) - kIp-cr], (A. 3) 
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C-sr = exp{iri(p +C L + P-cl + P+sl + P-sl) 
+ ni(-p +CR + P-cr -P+sr)], 

C +CL = esxp[ni(p +C L + P-cl + P+sl)], 

C-cl = exp[-7rip_cL], 

C+sl = exp[-7np_c L ], (A.4) 

C_ sl = exp [m ( -p+cL + P-cl - P+sl)}- 
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Figure Caption 



The effective potentials V as a function of 5a, for J" / J' = 0.3,0.6 and 0.9. It 
is observed that the minimum exists at 5a = at small values of J"/ J', whereas at 
some value of J" / J' the point 5a = becomes unstable. Phase transition occures 
and the system has an energy gap in the new phase. 
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